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In this paper we continue the analysis of quantum two-particle bound systems we have started in our earlier contribution 
(Kholmetskii, A.L., Missevitch, O.V. and Yarman, T. Phys. Scr., 82 (2010), 045301), where we re -postulated the Dirac 
equation for the bound electron in an external EM field based on the requirement of total momentum conservation, when 
its EM radiation is prohibited. It has been shown that the modified expression for the energy levels of hydrogenic atoms 
within such a pure bound field theory (PBFT) provides the same gross and fine structure of energy levels that had been 
predicted by the standard theory. Now we apply the PBFT to the analysis of hyperfine interactions and show the appear- 
ance of some important corrections to the energy levels (the 15-25 interval and hyperfine spin-spin splitting in positro- 
nium, 15 and 25-2P Lamb shift in hydrogen), which remedies considerably the discrepancy between theoretical predic- 
tions and experimental results. In particular, the corrected 15-25 interval and the spin-spin splitting in positronium prac- 
tically eliminate the available up to date deviation between theoretical and experimental data. The re -estimated classic 
25-2P Lamb shift as well as ground state Lamb shift in the hydrogen atom lead to the proton charge radius 
rp=0.841(6) fm (from 25-2P Lamb shift), and r p =0.844(24) fm (from IS Lamb shift), which perfectly agress with the 
latest estimation of proton size via the measurement of 2S-2P Lamb shift in muonic hydrogen, i.e. r p =0.84 184(67) fm. 
We also emphasize the universal character of PBFT, which is applicable to heavy atoms, too, and analyze 2S-2P interval 
in Li-like uranium. We show that the corrections we introduced provide a better correspondence between the calculated 
and experimental data than that furnished by the standard approach. The results obtained support our principal idea of 
the enhancement of the bound EM field in the absence of EM radiation for quantum bound systems. 
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1. Introduction 

The development of quantum electrodynamics (QED) in the second half of the past century and the 
unexampled success of this theory in physics of light hydrogenic atoms (where the deviation be- 
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tween theoretical predictions and experimental results gradually reaches relative values of 10" , 
and even 10" 13 ), completely convinced physicists about the outstanding correctness of QED with 
regards to the description of quantum phenomena. At the same time, modern progress both in the 
theory and experiment allowed recently the disclosure of a number of unexplained deviations be- 
tween theoretical predictions and experimental data in physics of light hydrogenic atoms, where the 
value of such deviation (A) substantially exceeds the corresponding uncertainty <j (both theoretical 
and experimental). Currently the revealed discrepancies between QED calculations and experimen- 
tal results are not the subject of wide scientific discussions, and thus it is worth to remind them be- 
low: 

- \S-2S interval in positronium (A/o»3.0) [1, 2]; 

- hyperfine interval in positronium (A/o»2.5) [1]; 

- proton charge radius r p derived from the classic 2S-2P Lamb shift and the ground state Lamb shift 
in hydrogen systematically exceed the value of r p obtained from particle physics (A/ cr varies from 3 
to 5 according to different estimations [3, 4]). 

One can add the result of very recent experiment [5], where the proton charge radius being 
estimated via the measurements of 2S-2P Lamb shift in muonic hydrogen (^=0.84184(67) fm), oc- 
curred to be less by 5.0 standard deviations than the modern CODATA value r p =0. 8768(69) fm [6]. 
Moreover, for muonic hydrogen the nuclear size effect contributes significantly (about 2 %) to the 
2S-2P Lamb shift and thus this new value of r p can pretend to be the most precise amongst all avail- 
able results. If so, its deviation from the value of r p extracted via the 2S-2P Lamb shift and 15* Lamb 
shift in the hydrogen occurs to be drastic. 

Thus, these facts may indicate on a possible presence of some still missed elements in the de- 
scription of quantum two-particle bound systems. 

In this respect we mention our recent paper [7], where we pointed out the known fact that 
quantum bound charges do not radiate at stationary energy states and thus their EM field consists of 
the bound (non-radiative) component only. This effect does not have a classical analogy, where, as 
known, an orbiting charge must inevitably radiate, and both bound and radiative EM field compo- 
nents equally participate in securing the total momentum conservation law. Hence the question 
emerges, which appears to be has not been asked before ref. [7]: How does Nature restore the ener- 
gy-momentum conservation law for quantum bound systems of charges, where their EM radiation is 
prohibited? 

As the general approach to the analysis of this problem, one can consider QED equations for 
Dirac field and EM field (the non-homogeneous wave equation for the four-potential [8]) and, using 
the formalism developed in ref. [9], to modify this wave equation for bound four-potential, cance- 
ling by such a way EM radiation. However, in the paper [7] we applied a much more simple but 
more illustrative way based on the analysis of one -body problem in the classical limit (still ignoring 
spin effects), via further taking into account the Bohr's correspondence principle. 

From the general viewpoint, it is clear that a possible way to implement the energy- 
momentum conservation law for an isolated system of charges with the prohibited EM radiation is to 
modify in an appropriate way their bound EM field (potentials) and/or the relationships between the 
fields (potentials) and their energy-momentum. Exploring this problem and considering the electron 
in an external EM field, we first have to distinguish its states with the total energy E>mc 2 (free elec- 
tron, where m is its rest mass and c the light velocity in vacuum), and E<mc (bound electron). 

At E>mc 2 the electron's energy spectrum is continuous, and its radiative EM field is not pro- 
hibited. In contrast, at E<mc the electron is characterized by the discrete energy spectrum, and it 
cannot emit EM radiation in a stationary energy state. Hence any possible modifications of the Dirac 
equation (DE) aimed to take into account the momentum conservation constraint for non-radiative 
nature of EM field of bound electron, should be made via the introduction of some step-wise func- 
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tion of the difference (E-mc ), which yields the common DE at E>mc . At the same time, it is rea- 
sonable to assume that for a bound electron (E<mc 2 ), the implementation of the total momentum 
conservation law in the absence of radiating EM field component would require the appropriate 
modifications of expressions for momentum of EM field and/or interaction EM energy. Thus the 
general form of DE for the electron in an external EM field, which reflects such modifications of 
bound EM field at E<mc 2 , reads [7]: 
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nonical momentum, e is the electron's charge, A b is the vector potential, <j is the Pauli matrix, <p is 
the scalar potential and the step-wise functions B n , r n are defined by the relationships 
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where b n , y n are some coefficients, we introduced. These coefficients are aimed to reflect the non- 
radiative nature of EM field of bound electron within the total momentum conservation constraint 
and thus, in the non-relativistic limit both b n and y n are put to be equal to unity. Hence their intro- 
duction in DE does not affect the gross structure of the electron's energy levels, characterized by the 
principal quantum number n. These coefficients, being constant for any fixed energy level, provide 
the Lorentz-invariance of the modified DE, when the non-relativistic limit is no longer assumed [7]. 
From the physical viewpoint, the presence of step-wise functions B n and r n in the Dirac equation 
reflects a discontinuity of the properties of the electron, which can emit EM radiation at the total 
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energy E>mc , but loses the ability to radiate at E<mc . 

In order to find the explicit form of the coefficients b n and y n for bound electron, we consi- 
dered in ref. [7] a model classical one -body problem, where EM radiation of the electron orbiting 
around a heavy nucleus, is prohibited, and further determined the way of modifying the bound EM 
field of this system, which guarantees the implementation of the total momentum conservation law. 
By such a way we found the classical limits for the coefficient b„ and y n , which occurred sufficient 
to find b n and y n themselves [7]. 

Deriving further the modified Dirac-Coulomb equation for the quantum one-body problem 
on the basis of modified DE, we arrived at the same gross, also fine structure of energy levels, as 
those furnished by the conventional approach, for hydrogenlike atoms, but obtained a small differ- 
ence in the value of average radius of electron's orbit [7]. A relative change of the radius of the elec- 
tron orbits for hydrogen-like atoms has the order of magnitude (Za) (where Z is the atomic number 
and a the fine structure constant), and is substantially smaller than the present experimental uncer- 
tainly in the measurement of atomic form-factors [1]. 

In the present paper we extend the approach of ref. [7] based on modified Dirac equation and 
named hereinafter as Pure Bound Field Theory (PBFT) to quantum two-body problems (section 2) 
and further analyze the hyperfine contributions to the atomic energy levels (section 3), which, in 
addition to the modified solution of quantum mechanical equations for the two-body problems, 
yields the replacement U^y n U'm the input of QED expressions for radiative corrections [7], without 
altering the core structure of this theory. We show that the corrections of PBFT to fine structure 
(having the order of magnitude of hyperfine interactions) are significant only for 1 S states of hydro- 
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genie atoms (sub-section 3.1). The correction brought by the PBFT to the hyperfine spin-spin inte- 
raction is completely negligible for the hydrogen-like atoms, except positronium, where such a cor- 
rection becomes significant (sub-section 3.2) to eliminate the available disagreement between calcu- 
lated and experimental data. In subsection 3.3 we derive the PBFT corrections to 15* and 25* Lamb 
shift, which substantially exceed the uncertainly (both theoretical and experimental) in their deter- 
mination for light hydrogenlike atoms. In section 4 we show that the corrections brought by the 
PBFT to common results practically eliminate totally the available up to date discrepancy between 
theoretical and experimental data for the \S-2S interval in positronium, spin-spin splitting in posi- 
tronium, classic Lamb shift and ground state Lamb shift in hydrogen. In particular, we derive the 
proton charge radius r p =0.837(8) fm (from 2S-2P Lamb shift), r p =0.840(24) fm (from IS Lamb 
shift), which completely agrees with the latest experimental result [5]. We also emphasize the uni- 
versal character of our approach, which is thus applicable to heavy atoms, too, and analyze the 2S- 
2P interval in Li-like uranium within the framework of PBFT (section 5), where the corrections we 
introduce provide a better correspondence between the calculated and experimental data than that 
furnished by the common approach. Finally, we conclude in section 6. 

2. Classical and quantum two-body problem 



In this section we consider first the hydrogenlike atom in the semi-classical treatment, suggested in 
[7], when the motion of electron with the charge e and mass m around the nucleus with the charge 
Ze and mass M is described in the classical way, but the EM radiation of this system is prohibited, 
like in the quantum case. For the one-body problem (M— »oo) we obtained the motional equation in 
the form: 
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which occurs sufficient to determine the classical limits of coefficients b„ (b = l + U/mc 2 ) and y n 

i 2 1 Ze 2 
( y - (l - v 2 / c 1 ) ) in the modified DE. Here U l = — \E e ■ E N dV = is the Coulomb interac- 
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tion EM energy for the one -body problem; E e and E N stand for the electric fields of electron and 
nucleus, correspondingly, and the integration is carried out over the entire three-space V. 

Analogously, when the mass of the nucleus M takes the finite value (two-body problem), we 
get the motional equation 

j t r ^ m+r M u l c2 ) v >» = -j t r ** ( M + r « u l cl Kf > ( 1 ) 

where 

U = ^\(E m -E M+ B m -B M )dV (2) 



v 



is the EM interaction energy. 

Eq. (1) allows defining the effective canonical momenta for the electron and nucleus 

Pbm =r m (m + r M U/c 2 )v m = y m m h v m = b m p m , p bM =y m (m + y m U / c 2 \ M = y M M h v M = b M p M (2a-b) 



where we introduce the effective rest masses 

m h = m(l + y M U/mc 2 ) = mb m , M h = M(l + y m U / ' Mc 2 )= Mb 



(3a-b) 



of both particles and designated p m ,pM the mechanical momenta of electron and nucleus, correspon- 
dingly. Herein, like in ref. [7], we supply the quantities of PBFT by the subscript "6", and introduce 
the factors 



b = 



it 2 

^ mc J 



l+ /m 



v 



Mc 1 



2 / 2 



(4a-d) 



Then the Hamilton function, written in the weak relativistic limit to the accuracy c , becomes 
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The approach consisting in going from the classical description to quantum description of 
two-particle bound system and taking into account the Dirac equation, is not directly applicable to 
the two-particle case, where we should address either to the Bethe-Salpeter equation, or the Breit 
equation without external field [10], or to their appropriate modifications. Though the Breit equation 
is not fully Lorentz-invariant and represents an approximate, it is the most convenient and illustra- 
tive for the analysis of PBFT corrections, resulting due to the replacements 



m^b mn m, M^b Mn M, U^y mn y M nU, 



(5a-c) 



where the classical limits of the introduced coefficients b mn , b Mn , Ymn and y Mn are determined by the 
respective equations (4a-d). Besides, we define the operators of canonical momenta 
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Pbm = PbM = Pb = -iKV , where we take into account that due to the total momentum conservation 
law, Pbm = -PbM =Pb- Hence the Breit equation modified by the replacements (5a-c) reads: 

(H b0 +H hl +... + H b5 )%{r m ,r M )= W{(r m ,r M ), (6) 



where %{r m ,r M ) is the wave function having 16 spinor components, r m , r M are the position vectors 
for each particle, W is the energy, and the Hamiltonian components H b i have the form: 
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where * = -= — , and r=r m -r M . 



Here the first operator H b0 represents the counterpart of PBFT to the conventional 
Schrodinger operator; H bl is the relativistic expansion of H b0 ; the operator H b2 takes into account 
the retardation in the interaction of electron and nucleus; H b3 describes spin-orbit interaction; H bA 
is responsible for the contact interaction, and H b5 stands for spin-spin interaction, where we take 
into account that in PBFT, the interaction of two spins (currents) is y mn y Mn times larger than in the 

common theory, due to the modified field equation for magnetic field and Lorentz force law in pure 
bound field electrodynamics (eqs. (13d) and (11) of ref. [7], correspondingly). The same factor 
YmnYun i s introduced into the term of spin-orbit interaction, because the spin-orbit term originally 
includes the electric field of the nucleus Eh, which is related to the usual electric field E=Zerlr by 
the relationship E h -^y mn y Mn E, implied by eq. (5c). (In fact, the substitution E b -^y mn y Mn E originates 
from the field equations and Lorentz force law applied to two-body problem in pure bound field 
electrodynamics (see eqs. (13a) and (11) of ref. [7]). 

For further analysis it is convenient to reduce the Breit equation to the Schrodinger-like type 
by analogy with ref. [11]: 
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In order to solve eq. (7), it is convenient to apply the substitution 

r = r'/{b mn b Mn y mn y Mn ), (9) 



which will allows us to present the Hamiltonian in eq. (7) as the sum of non-relativistic Schrodinger 
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term and perturbation. Indeed, taking into account that p h = —h V \ = —b mn b Mn y mn y Mn fi V r , , we 
transform eq. (7) as follows: 
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The obtained eq. (10) completed by the expressions (8), (9), (11) represents the basic equa- 
tion for the quantum two-body problem within the framework of PBFT. Here one should recall that 
eq. (10) itself, like the source Breit equation, is semi-relativistic, and it is valid to the order (Za) 4 . At 
the same time, the factors b mm b Mn , Ymn and y Mn , being explicitly determined to the orders {Za) and 
(Za) 4 (see below), allow us to analyze the specific PBFT corrections to the order (Zaf, which cor- 
responds to the scale of hyperfine interactions. The determination of these corrections is the next 
goal of our analysis, but, first of all, let us show that eq. (10) yields the same gross and fine structure 
of the atomic energy levels, as the one furnished by the common approach. 

In the zeroth approximation, when the terms of order (vie) and higher are ignored, we get 
from eq. (10) the Schrodinger equation expressed in r' -coordinates: 
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where m R = mM /(m + M) is the reduced mass. Hence we obtain the well-known solution 

m R c 2 {Zaf 



W n = 



2n 2 



along with the common Schrodinger wave function expressed via r' -coordinates. This result allows 
us to obtain the coefficients b mn , b Mn , Ymn, YMn at least to the order (Za) , based on their respective 
classical limits (4a-d), taking into account the known relationships (e.g. [12]): 
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where vr is the reduced velocity. Hence, via the comparison of eqs. (12a-d) with eqs. (4a-d), we ob- 
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tain the factors b mn , b Mn , Ymn, YMn to the accuracy (Za) as follows: 
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Further on, using eqs. (13a-d), we derive the product 
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to the accuracy of calculations (Zd) . 

Applying equations (13), (14), as well as the equality [11] 
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Substituting this equality into eq. (10), and ignoring the PBFT factors b m „, bMn, Ymn, Yun in the terms 
of the order (Za) 4 , we obtain: 
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where the term U[p m ,p M ,r') differs from the term (8) by the omission of PBFT factors b mn , b Mn , 
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Excluding further the term of spin-spin interaction in the expression for L/(p m ,p M ,r') (last 
term in the rhs of eq. (8)), we arrive at the common solution for the Dirac-Recoil (DR) contribution 
to the energy levels, written to the order (Za) 4 [3]: 
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and j is the quantum number of total angular momentum (j=l+s, I is the angular momentum, and s 
the electron's spin). 

Thus the corrections of BPFT to the common solutions of equations of atomic physics may 
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emerge at least in the order (Za) 6 , which corresponds to the scale of hyperfine interactions, and 
which will be determined in the next section. As we will see below, these corrections considerably 
eliminate the available discrepancy between theory and experiment in physics of light hydrogenlike 
atoms. 



3. Hyperfine contributions to the atomic energy levels: Pure 
Bound Field theoretical approach 



We have shown above that the solution of the modified Breit equation gives the same gross as well 
as fine structure of energy levels for hydrogenlike atoms, as that yield by the conventional solution, 
up to the order of (Za) 4 . In this section we analyze hyperfine contributions to the atomic energy le- 
vels in PBFT, which can be presented in the form 

(W b ) n y = {K R ij + (W b ) HFS + (kl, , d9) 

where \W b ) HFS is the hyperfine splitting due to spin-spin interaction, and (L b ) nlj is the Lamb shift. 

Herein the subscript "Z>" reminds that the corresponding energy term is evaluated within the PBFT 
framework. 

In what follows, we consequently analyze the fine structure corrections (sub-section 3.1), 
corrections to spin-spin interaction (sub-section 3.2), and corrections to the Lamb shift in light hy- 
drogenlike atoms (sub-section 3.3) within the approach of PBFT. 

3.1. Corrections to the fine structure 



In this sub-section we show that the corrections of PBFT to fine structure of light hydrogenlike 

atoms with m«Mhave the order of magnitude of [Zaf — and higher, and due to their scaling as 

M 

n' 6 (or n' 5 ), they should be taken into account practically only for the ground states. In the case of 
positronium, the fine structure corrections occur significant not only due to the equality m=M, but 
also due to the PBFT correction to the annihilation term. 

First of all, we extract from the Hamiltonian of eq. (10) the terms of fine interactions with the 
order (Za) 4 properly modified in PBFT, omitting at this stage the contribution due to spin-spin inte- 
raction. Hence, via eqs. (8) and (10) we obtain the operator of fine interaction in the form 
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is the term of contact interaction, and 
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is the term of spin-orbit interaction. 

We point out that the operators (21)-(23) are presented in r-coordinates, whereas in the Ha- 
miltonian (10) they should be expressed through r'— coordinates. In the latter case, each of the terms 
(21)-(23) can be calculated with taking into account of eq. (9), as well as the equalities 



Phi 1 "') = Knb M nY m nYMnPt( r ) > K r ) = Kn h Mn Y mn Y Mn 8 ^f) ■ 
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Introducing the replacements (24-25) into eq. (21-23), and expressing the terms (22), (23) via the 
fine structure constant, we obtain the operator of fine interactions in r' -coordinates as follows: 
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Substituting this value in eq. (10), we get 



n 2 vlb^ n z v z M 



2w2 ; 



r Ma 



Ze 



2m 



2M 



fine 



¥ {r') = 



W 



h h 22 

inn MnYmn Y Mn 



v(4 



(27) 



We emphasize that to the order (Za) , eq. (27) is equivalent to eq. (17), which yields the 
same fine structure of energy levels, like the common approach. Now our goal is to determine the 
fine structure corrections of PBFT in eq. (27), which may emerge in the order (Za) 6 . To solve this 
problem, we need to determine factors b mn , bMn, Ymn, YMn to the order (Za) 4 . 

For this purpose, first of all, we find the classical coefficients b m , b M , and y m , y M to the accu- 
racy c' 4 based on their definitions (4a-d). Concerning the factors b m , b M , their modification issues 
from the two circumstances: 

- the factors y m and y M in eqs. (4a-b) are no longer adopted to be equal to unity, but they have to be 
determined to the accuracy (vie) ; 

- in the expression for interaction EM energy (2), the term of magnetic interaction energy is no 
longer ignored. Taking into account that for the bound EM field B = (vxE)/c, we obtain for the 
circular motion in the classical two-body problem 

R R _ D D _ V m V M p p mM V R p F 

D m' D M- D m D M- 2 m M / , ,\2 2 m ' M ' 

c (m + M) c 



Substituting this equality into eq. (2) and following the quantum mechanical definition of U , we 
obtain 
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1 + 



Mn 

mc z 



(Zaf ( 1 | (Za) 



m 



2n 2 (m + MfX (m+M) 



1 + 



mM (Za) 



2\ 



M 



(m + M) 



M 



(Zaf m 2 M (Zaf mM 2 



n 2 m + M 2n A (m + M) 3 



(m + Mf 



(28a) 



to the accuracy (Za) , where we have used eq. (13d), putting 



\_(Zaf 



m 



(m+M) 2 



-1/2 



m 



2n (m + M) 2 



with the sufficient accuracy of calculations. 
Similarly we determine the factor 



_^ (Zaf m (Za) 4 M 2 m (Za) 4 



Mm 



n 2 m + M 2n 4 (m + M) 



(m + M) 3 



(28b) 



We point out that expanding the classical coefficients b mn and 6m« (eqs. (4a-b)) to the accura- 
cy (ylcf, we adopt that the involvement of non-Coulomb interactions does not affect their values. 
Indeed, at the semi-classical level, as shown in ref. [13], such non-Coulomb interactions are exactly 
counteracted by corresponding change of the Coulomb interaction energy due to the proper variation 
of the radius of electron's orbit. Hence the resultant action of Coulomb and non-Coulomb interac- 
tions is equivalent to the Coulomb interaction alone with the fixed electron's orbit. This observation 
concurrently implies that the overall change of the energy of semi-classical system "electron plus 
nucleus", for example, due to spin-orbit interaction exhibits as a proper change of kinetic energy of 
the orbiting electron and nucleus by the energy of non-Coulomb interaction [13]. Therefore, the tak- 
ing into account of fine interactions does modify the factors y mn y Mn in the orders (Za) 4 and higher. 
In order to find the related corrections, it is convenient to use the relationships between the Lorentz 
factors and momenta of particles 



Yrr, 



1 



l-p 2 b /m 2 c 2 



7a 



1 



\-pIIm 2 c 2 



where pt is the value of non-relativistic momentum of electron (nucleus), and to the order (Za) the 
averaged value of p 2 h is defined by eq. (15). 

With the involvement of fine interactions, this expression is modified as 



p 2 = 2m { 



{w+v fne )+ 



Ze 2 



Ze 2 



Taking into account that W H = mc 2 (Za) 2 [1 1], we obtain 
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Y mn 



1 



1 



1 (Zaf M 2 2m R (V h (r%,/ /M » " ( Za ) 2 m 2 2m R (V h (r%, 
n (m + Mf mc n (m + Mf M 2 c 



.(28c-d) 



Having obtained the coefficients b mn and bun (eqs. (28a-b)) and y mn , yu n (eqs. (28c-d) we now 
in the position to calculate the fine structure corrections of PBFT to the order (Za) 6 on the basis of 
eq. (27). In particular, we find that 



P%n , P\ 



■ + - 



2m 2M 2m 



P 2 _ m „ c ^^^^ n(r .) L %I^L 



R 



n 



(m + M) 2 



n 



(m+M) 2 



(29) 



Further, we determine the product 



W 



m n Mn Y mn Y Mn 



= w- 



m 



W 

2 2 ' which for the nS state is equal to 

^ mtP MnY mn Y Mn 

, 6 .„2 i y-2 



R c 2 (Za) 4 Mm — (Zaf m R c 2 (Za) 6 5mM m R c 2 (Zaf m 2 M 



(m + M) 



- + V f 



n 2 An (m + M) 2 



(30) 



Here we have used eqs. (28a-d), and in the terms of the order (Za) we put W=W( 



2n(m+M) 4 

m R c 2 (Za) 2 



0«" 



2n 2 



2 — 

while in the terms of the lower order (Za) we presented W = W Qn +Vf me . 

In order to find the PBFT corrections for the term (V h (r } )) fine , we can write it in the form 

{vXr% ine ={v(r% ine +5V fine , 

where (v(r')) fme represents the common operator of fine structure, and 8V fme contains the specific 

corrections of PBFT determined by the substitution of factors b mn , b Mn , Ymn, Ymh in eq. (26). Since the 
operator (v(r')) fine itself has the order (Za) 4 , it is enough to use the PBFT factors (13a-d) written to 

the accuracy (Za) . In particular, for nS-state of hydrogenlike atom, the straightforward calculations 
yield: 



8V r , 



M 2 +m 2 (Zaf -- | 2mc 2 (Zaf 
(M+m) 2 n 2 n(M+mf 



(mM 4 -m 2 M 3 +m 3 M 2 )- 



9 mc(Za) 6 mM 2 (M 2 +m 2 ) 



8 n 6 



(M + mf 



(31) 



Substituting eqs. (28-31) into eq. (27), we obtain after lengthy, but straightforward calcula- 
tions for nS-states: 



2m R r 



W 



2mc 2 (Za) 6 ( 2mM * + m i M ^ +2 m'M 2 )+ 



mc 2 (Za) 6 mM 2 + m R c 2 (Zaf m 2 M 2 



n~(M + m) 5X ~'" ' '" " n 6 ^ m + M J 2n 6 (m + M) 

2 (Zaf (M 2 +m 2 )mM 2 9 mc 2 (Za) 6 mM 2 (M 2 +m 2 ) 
+ THC 2n 6 (m + M) 5 8 n 6 (M + m) 5 



(32) 



y,(r') 
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Designating the term in the bracket of rhs of eq. (32) as W common , we rewrite this equation in the 



short form 

2 v„2 



p' Ze 



2m R 



which represents the common Breit equation without spin-spin term written in the Schrodinger-like 
form (see, e.g. [11]). Therefore, the PBFT correction to the atomic nS state energy levels is equal to 

5W PBFT = W - W_ = 2 ™f} Za l {^M 4 + m 2 M> + 2m' M 2 )- 

n(M + m) 

mc 2 {Zaf mM 2 m R c 2 {Zaf m 2 M 2 2 {Zaf (m 2 +m 2 }nM 2 9 mc 2 (Za) 6 mM 2 (M 2 +m 2 ) 
n 4(m+M) 3 2n 6 (m + M) 4 2n (m+M) 5 8 n (M + mf 

(33) 

One can see that for atoms with m«M, the correction (33) has the order of magnitude 
mc 2 (ZaJ > m/M and scales as n' 6 or n' 5 . Therefore, there is no practical need to calculate this correc- 
tion for n>2 and /^O. 

In fact, the correction (33) is significant for IS state only, and for the hydrogen atom it is 
equal to 

l imc 2 {Zaf £ =14.8 kHz (34a) 

Note that this correction is positive, and it reduce the value of 15* state (which is negative). For 2S 
state of hydrogen we obtain from eq. (31) SlV^ ne \2S) =0.91 kHz. Hence the correction to lS^S in- 
terval is equal to 

StV^ ne (lS-2S) =13.9 kHz. (34b) 

Below the correction (34b) will be involved into the re-estimation of 15* Lamb shift in hydrogen. 

For 15* state of muonium, this correction is about 0.13 MHz, and is few times less than the 
theoretical uncertainty in calculation of \S-2S transition (-0.30 MHz [3]), and much less than the 
experimental uncertainty in the measurement of this interval (9.8 MHz [14]). 

For 15* state of positronium (m=M), the correction (33) becomes 

W;;SlS)= 2l ™^ a)6 =3.07 MHz, 

and for the \S-2S interval, as follows from eq. (33), 

5Wj* e (IS - 2S) =2.95 MHz, (34c) 



which will be involved below into the re-estimation of 1S-2S interval in positronium. 

Besides, for positronium the Breit potential includes the additional annihilation part (e.g., 
[10, 11]), which in PBFT acquires the form 
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b y 

n t n 



^|^(3+ ° + .a_)s{r) 
y 2m b n c 



where a + ( a ) belongs to positron (electron), and we designated b mn =bMn = b n , Ymn = YMn=Yn- In r } - 
coordinates, due to eq. (25), this operator acquires the form 

V 2m c 



We average this term with the wave function for /=0 [11] 



Ho) =4 



4n 



' Zam R A 
{ n 2 J 



3/2 



and taking into account that for positronium b n = 
for m=M), we derive in the orthopositronium case 

w = 

V o J ann 



\ (Za) 
2n 



2\ 



' y mn 



\_{Zaf 
4n 2 



-1/2 



(see eqs. (13) 



■-<x) u 



mc 2 (Za) 4 j 3(Z«) 



2^ 



j j 



Hence the correction of PBFT to annihilation term reads: 

3mc 2 (Zaf 



5W n 



16w- 



(35) 



The correction (35) decreases the value of \S-2S interval in positronium by 3.40 MHz, and thus, it 
should be added to the fine structure correction (34c). A comparison of the theory and experiment 
for \S-2S interval in positronium will be done below in sub-section 4.1. 

3.2. Corrections to hyperfine splitting of energy levels due to spin-spin interaction 

Now we analyze the contribution of spin-spin interaction into the Breit potential, which in PBFT has 
the form 



fcML = 



KnKnYmnYMn ^b^Mb^S 



z e 2 h 2 y mn y M n ( %r % 3 (<yX%- r ) 8;r j( r ) 



(36) 



(the last term of eq. (8)). Being expressed via r' -coordinates, this operator reads: 



( V b( r ')l- S =b mn b M J m ,!r j 



e 2 K 



/! AmMc 2 



° m - °m a (q„- *•')(% -r') 8tt 



j3 



S{r<) 



(37) 



where we have used eqs. (9) and (25). Designating 
(V(r')\ 



AmMc 2 



e 2 h 2 ( «_ 3 MM_fe^ 

3 U 



v r r 

(the common Hamiltonian of spin-spin interaction expressed via r' -coordinates), and substituting 
expressions (13a-d) for PBFT factors, determined with the sufficient accuracy {Za) , we obtain 
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n 2 (M + mf 



(38a) 



This relationship is also valid for the energy of spin-spin interaction, obtained via the averag- 
ing of operators (V b (r')) s _ s and (v(r')\_ s with the wave function y/(r'f. 



1 (Zaf 2mM 
n 2 (M + mf 



(38b) 



Now it is important to remind that the energy W s . s contains the ratios of magnetic moment to 
mass both for the electron and the nucleus, which are determined experimentally by means of the 
Zeeman effect. Since in PBFT the operator of interaction of magnetic dipole with an external mag- 
netic field is, in general, modified, the appropriate corrections to the measured values "magnetic 
moment/mass" ratio should be clarified, too. 

As known, the operator of interaction of two bound particles (electron and nucleus) with the 
external magnetic field reads [1]: 



Knag Sm i S m ' Sm ^ M ' ^' ' 



(39) 



where g m , g M are the g-factors for bound electron and nucleus, correspondingly. Being added to the 
Breit operator of eq. (10) with the PBFT corrections m^b mn m, M—>b mn M, B^/ mn / M nB, this operator 
acquires the form: 



(n) = 

V b //naff 



1 



2 2 



^mn^MnY mn Y Mn 



eh / s_ Zeh / x 

Sm ^ , Y mnY Mn\ S in ' " ) Sm , „ Y mnY Mn\ S M ' "I 



2mb 



2Mb 



Mn 



Averaging this operator with the Schrodinger wave-function yKf), due to the normalization require- 
ment y/(r) = {b mn b M j m j Mn y l2 y/{r') implied by eq. (9), we obtain 



r \ 

V J 



/ jrr \ 1 1 2 z 

\ b /mag MnY mn Y Mn 



mag 



S, 



2m 



2M 



(40) 



where (W h ) mag gives the Zeeman splitting of energy levels in PBFT framework. Herein in averaging 



of 



/ A 

V t 
\ J 



we put B(r) = const within the atomic scale, which is always fulfilled in practice. 



mag 



Inserting eqs. (13a-d) into eq. (40) and putting Z=l (which is enough for our immediate pur- 
poses), we obtain for nS-state 



(w\ f i ( Za ) 2 2Mm 1 w ( Za Y efl (( - \ r\ 1 

WbLag- 1 n i ( m+M y w »™g „2 ^y\s m s m Sm*m )■ * ) A/f ^ 



2 ' ' mag 
) 



n 2m 



M + m 



(41) 



where W mag stands for the Zeeman splitting of energy levels, obtained via the averaging of common 
operator (39). For the sublevel F=\ used for the measurement of Zeeman effect, s m - s M =0, and eq. 
(41) acquires the form 
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M 



2m 



1 



M + m 



(42) 



where we supply the mass M by the subscript "Z" ("Zeeman effect"), in order to distinguish it from 
the mass M in eqs. (38a-b), designating the mass of the nucleus in the measurement of spin-spin 
splitting. Thus, the ratio of "magnetic moment/mass" derived from the Zeeman splitting should be 
corrected in PBFT for each bound particle with taking into account of the relationship (42) between 

For the 15* state of hydrogen and heavier atoms, the analysis of this corrections in the estima- 
tion of spin-spin interval (38b) is not practically important, because the term [Zaf m ^ W s _ s of 

(M + m) 

eq. (38b) itself is less than 100 Hz and is many times smaller than the nuclear-structure corrections 
to the 15* hyperfine splitting, which vary from tens to hundreds kHz [15-17]. Thus for hydrogen the 
PBFT correction to spin-spin splitting occurs negligible, and we put 



(43) 



within the range of the present uncertainty in calculation of Wf_ s . 

Considering the spin-spin splitting of 15* state of muonium, we can also ignore the correction 
to the ratio "magnetic moment/mass" for the electron, because one can show that it induces the 
PBFT correction of the order 100 Hz, which is much smaller than the present theoretical uncertainly 
in calculation of spin-spin interval in muonium (about 500 Hz [1]). Further, for muonium we can put 
with a high accuracy g m =gM in the term containing (Zd) . We also use the known fact that the ratio 
"magnetic moment/mass" for bound muon is determined with the best accuracy via the Zeeman ef- 
fect in muonium [1], so Mz=M. With these equalities eq. (42) yields 



( W y* = 

V o /mag 



\_{Zaf 



2mM 



(m + Mf 



W, 



Mu 

mag ' 



(44) 



which shows that the PBFT correction to spin-spin interval (38b) has exactly the same structure as 
PBFT correction to the "magnetic moment/mass" ratio for bound muon (44). As a result, both cor- 
rections exactly compensate each other, and we get 



1- 



(Zaf 2mM 
n 2 (M + m) 2 

(Zaf 2mM z 
n 2 {m + M z f 



(45) 



For positronium we have quite different situation (M=m, but Mz~ 10 m if the magnetic mo- 
ment/mass ratio for electron is taken from the Zeeman splitting in hydrogen or heavier atoms), and 
the correction to spin-spin interaction (38b) dominates over the correction to the "magnetic mo- 
ment/mass" ratio for bound electron. Hence we involve the correction of eq. (38b) solely, which 
yields 
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This equation will be used in sub-section 4.2 for PBFT correction of spin-spin interval for IS 
state of positronium and for comparison with modern experimental data. 

3.3. Corrections to the Lamb shift. 

The corrections of PBFT obtained above in sub-sections 3.1 and 3.2 originate from the appropriate 
modification of DC equation suggested in ref. [7]. Analyzing now radiative corrections to the atomic 
energy levels, one should emphasize that the modifications in the Dirac equation do not touch the 
core structure of QED, but at the same time imply the replacements (5a-c) in the input of QED ex- 
pressions. On the basis of this result we derive below the corrections to the Lamb shift L for light 
hydrogenlike atoms, which emerge in PBFT. 

It is known that the dominant terms of the Lamb shift arise due to a finite radius of the elec- 
tron (r 2 ^ , which continuously emits and absorbs virtual photons, as well as due to vacuum polariza- 
tion. 

The finite radius of the electron induces a deviation from the Coulomb potential [3] 

= Ur 2 )AU « f h(Z«r ™ = -MZa)- lZ ^5{r), (47) 
6 x ' in mi m 

where A is the Laplacian. According to eqs. (5a-c), for the bound electron the mass m is replaced by 
b mn m, while the Coulomb potential U is replaced by U'=y mn y Mn U~y mn U, where we can put ^mh=1 
with the sufficient accuracy of calculations; further we also put 6m«=1 in the correction to the Lamb 
shift. Thus the latter equation acquires the form 



SV'r^ radlU s= rfyln(Za)- 2 -^d{r) = f^ fmite radlus , (48a) 
3Z> m h. 



where SV {miie radius is defined by eq. (47). 

The contribution due to vacuum polarization [3] is also proportional to AU/m , and thus the 
latter equation remains in force for this correction: 

'^polarization — ~ 2~ '^'polarization ■ (48b) 

The total contribution SV Xo ^ is defined as the sum SV M = SV fimte radms + <^ polarizatlon , so that for the 
total perturbation we get 

^total= fractal • (48C) 
Kn 

The correction to the energy level is given by the matrix element of the total perturbation (48c), 
where we need to take into account that due to normalization requirement implied by the scaling 

transformation (9), y/(r) = {b mn b Mn y mn y Mn f 2 y/(r') ~ bj 2 yj 2 y/{r' ,3,<p). Hence 
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AW PBFT = {nS\SV'JnS) = b n J r , 



3 3 



= b y 4 AW = y 2 AW 

u mnl mn I mn L - y " ■ 



where AE denotes the value of energy shift, obtained within QED, and we have used the equality 
b mn Ymn = 1 ' followed from eq. (14) at M— >co. 

Thus the Lamb shift at the given energy level corrected within PBFT, reads: 



(O, 



nlj ) PBFT ' mn n H 



where L nlj stands for the Lamb shift calculated in QED. 

For the 2Sm-2P\a Lamb shift both levels have a principal quantum number n=2, and we get 

( L b ) 2S -2P = Ym2 L 2S-2P = L 2S-2P (l " ( Z «)' / 4 ) ' • (50) 



Thus the correction induced by PBFR to the 2S-2P Lamb shift is equal to 

3L2S-2P = ih \s-2p ~ L2S-2P = L 2 s-2p (l ~~ {Zoc) I A j — 1 



(50a) 



Numerically this value is equal to 13.8 kHz, which exceeds substantially the measured precision for 
Doppler-free two-photon laser spectroscopy [18]. 

Below we compare the experimental and theoretical values for the corrected 2S-2P Lamb 
shift in hydrogen (sub-section 4.3) and 2S-2P Lamb shift in He + (sub-section 4.4). 

Eq. (49) is, in general, also applicable to the 15* Lamb shift Lis in hydrogenlike atoms. How- 
ever, its direct measurement is impractical until effects of nuclear structure are known accurately 
enough. In order to eliminate the influence of these effects, the data at least of two measurements are 
involved: for hyperfine intervals in the ground state and metastable states (for example, for the 1 S 
and 25* states). Since the bulk contribution to the Lamb shift scales like n , then the difference 
8(^/) 2 5 ~^Wh P f)\s allows us canceling substantially various contributions caused by the short dis- 
tance effects. However, the factors y mn differ from each other for 25* and 15* states, and calculation of 
the corrected 15* Lamb shift {L b \ s is not straightforward. 

In order to introduce the PBFT corrections to the IS Lamb shift, to be convenient for practic- 
al applications, one need to look closer at the typical methods for its theoretical estimation. In prin- 
ciple, the 15* Lamb shift could be extracted from the experimental data on the transition frequencies 
between the energy levels with different numbers n. One should emphasize that the intervals of 
gross structure are mainly determined by the Rydberg constant R. In order to disentangle measure- 
ment of the 15* Lamb shift from the measurement of the Rydberg constant, one can use the experi- 
mental data on two different intervals \S-2S and 2Si/2-SD 5 / 2 of hydrogen [19]. Theoretically these 
intervals can be presented as [3] 

E\s-2s — (^25 V2 _ ^i5 l/2 )+ As V2 _ As V2 ' (51a) 

^2S-SD — ty&D y2 ~ WlS l/2 )+ L&D 5/1 ~ ^2S l/2 ' (51b) 

where W% R is the leading Dirac and recoil contribution to the position of the respective energy level 
(eq.(18)). J 
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The differences of the leading Dirac and recoil contribution in the rhs of eq. (51) are pro- 
portional to the Rydberg constant R plus corrections of order a R and higher. One can construct a 
linear combination of these intervals which is proportional to a 2 R plus higher order terms 
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(52) 



Then the difference of the leading Dirac recoil contribution in the rhs of eq. (52) can be calculated 
with a high accuracy, due to the suppression factor a 2 , and it practically does not depend on the ex- 
act value of R. Hence the linear combination of the Lamb shifts in the rhs of eq. (52) does not de- 
pend on R, too. The bulk contribution to the Lamb shift scales as \ln which allows using the theo- 
retical value L SDs2 =71.51 kHz [3] without loss of accuracy. The 2S Lamb shift can be extracted 

from the data on the classic 2S-2P Lamb shift, so that 
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(53) 



Herein in rhs, the first term in square brackets is computed theoretically, the second term in square 
brackets is determined experimentally, while the last term is extracted from the data on the classic 
2S-2P Lamb shift. Within PBFT, the first computed term in the rhs of eq. (53) should be corrected 
by adding the fine structure correction (34b). Besides, one has to correct within PBFT the 25* Lamb 
shift, using the data on the classic 2S-2P Lamb shift. Hence the expression for the 15* Lamb shift for 
the hydrogen within PBFT acquires the form 

[L b ) ls [kHz] = L ls [kHz]+ 13.9 + jSL 2S [kBz] ; (54) 

the PBFT correction to 25* Lamb shift SL 2S will be found below in sub-section 4.5, where the value 
(54) will be calculated. 

4. Corrections to the atomic energy levels and comparison with 
experiment 

In this section we analyze the hyperfine contributions to the energy levels of light hydrogenlike 
atoms, where the discrepancy between theory and experiment exceeds the uncertainties in their de- 
termination and apply the appropriate corrections of PFBT derived above. We show that the correc- 
tions of PBFT provide a perfect conformity between theoretical and experimental values for all pa- 
rameters listed in the introduction section: \S-2S interval in positronium (sub-section 4.1); spin-spin 
splitting in positronium (sub-section 4.2), proton charge radius derived from the classic 2S-2P Lamb 
shift (sub-section 4.3), proton charge radius derived from the ground state Lamb shift in hydrogen 
(sub-section 4.5). We also pay a separate attention to the 2S-2P Lamb shift in He + (sub-section 4.4). 

4.1. 1S-2S interval in positronium. 

Modern theoretical value of this interval is [1] 

£, P /_ 2S = 1 233 607 222.2(6) MHz, (55) 
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and the most precise result of experimental measurements is as follows: 

1 233 607 216(1) MHz [2]. (56) 

One can see that the deviation between the values (55) and (56) more than three times larger than 
the uncertainty of measurement of 15-25' interval. 

Now we introduce the PBFT correction to \S-2S transition as the sum of eqs. (34c) and (35): 

(ml:^S-2S) = SfV h P %lS-2S) + (3V h ) ann =6.35 MHz. 

Hence the 1S-2S interval in positronium corrected in PBFT becomes 

{E b J;_ 2S =E 1S _ 2S -(A l$ s s _ 2 =l 233 607 215.8(6) MHz, (57) 

which perfectly agrees with the experimental value (56). 

4.2. Spin-spin interval in hydrogenlike atoms 

In sub-section 3.2 we have shown that the correction of PBFT to hyperfine spin-spin interaction 
occurs quite negligible for atoms with m«M (e.g., hydrogen, eq. (44) and muonium, eq. (45)). For 
positronium we derived eq. (46), which now will be used for the comparison with experimental data. 
The theoretical value of hyperfine splitting in positronium is [1] 

Wf_\ =203 391 .7(8) MHz, (58) 

which does exceed the corresponding experimental data 203 389(2) [20] and 203 387(2) [21]. 

Eq. (46) allows us to compute the corrected PBFT value of hyperfine spin-spin interval in 
positronium, using the numerical value (58): 

(W b f_ s =203 386(1) MHz. (59) 

This result is already in a good agreement with the experimental data. 

4.3. 2Si/2-2P 1/2 Lamb shift in hydrogen. 

It is well known that the dominant problem of exact theoretical evaluation of the classic Lamb shift 
L 2S _ 2P is the uncertainty arising from the proton charge radius r p . Due to this reason many authors 

reverse the problem, and estimate r p from the obtained data on L 2S1P shift (see, e.g., [1]). It is also 

known that the estimated value of r p via the measurement of classic Lamb shift systematically ex- 
ceeds the magnitudes of r p , obtained in the electron-proton scattering data and other methods for 
evaluation of r p in physics of elementary particles [4]. This prompted scientists to assume [22] that 
the uncertainties in estimation of r p in the experimental particle physics are significantly underesti- 
mated. However, the very recent estimation of proton charge radius via the measurement of 2S-2P 
Lamb shift in muonic hydrogen gives the value r p =0.84 184(67) fm [5], which is substantially lower 
than the CODATA value r p =0. 8768(69) fm [6]. It is also important that for muonic hydrogen the 
nuclear size effect contributes significantly (about 2 %) to the 2S-2P Lamb shift and thus this new 
value of r p can pretend to be the most precise result amongst all published. 
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Below we will show that the PBFT correction (50a) to the 2S-2P Lamb shift removes the 
exiting remarkable disagreement between the estimation of r p from the classic Lamb shift data and 
estimation for muonic hydrogen. 

First we determine the factor y m 2, which for hydrogen and muonic hydrogen atoms has the 

value y m2 - (l-a 2 /4) 1/2 =1.0000066. For muonic hydrogen, where the nuclear size effect contri- 
butes significantly to the total 2S-2P energy interval, the corrected Lamb shift (50) with the factor 
Ymi computed right above does not practically affect the proton charge radius estimated in ref. [5]. In 
particular, using the parameterization (1) of ref. [5] for the 2S-2P energy difference, one can show 
that the correction (50a) influences the estimated proton size in the order of magnitude 10" 4 fm, 
which is below of the measurement uncertainty [5]. In contrast, for the hydrogen atom the correction 
(50a) and finite nuclear size effect have comparable values, and the proton charge radius derived 
with and without correction (50a) acquires a difference to be substantially larger than the meas- 
ured/calculated uncertainty. 

In order to estimate the proton charge radius from the classic Lamb shift, we use the parame- 
terization 

L 2S _ 2F (r p )=A + Br p 2 , (60) 

which is based on the known fact [3] that the term proportional to r p is additive. Here A and B are 

the coefficients, whose numerical values can be found via common calculation of 2S-2P Lamb shift 
in hydrogen [3] for different values of proton charge radius [23, 24]: 

^=1057695.05 kHz, 5=195.750 kHz/fin. (61a-b) 
In the framework of PBFT, the eq. (60) is appropriately modified: 

( L 2S-2p)pbFT = Ym2 A + f^D/ » ( 62 ) 

where (r PBFT ) is the proton size predicted by PBFT. Putting in eq. (60) 

r p =0.876(6) fm (63) 

(COD ATA value [6]), and equating (60) and (62), we derive the expression for the new proton size 
evaluated in PBFT: 




with the sufficient accuracy of calculations. Using the numerical values (61) and (63), we obtain 
(w), = 0-834(6) fm. 

This estimation is much closer to the proton size derived in ref. [5], than the CODATA value 
[6]. At the same time, now we recall that the modern CODATA value (63) incorporates the experi- 
mental data in both particle physics and atomic physics, and, in general, is less than the proton size 
derived from the classic Lamb shift solely. In particular, the modern data on 2S-2P Lamb shift in 
hydrogen obtained by various authors within the common approach (see refs. [1,3] and references 
therein) define the range of variation of the values of r p between 0.875 fm and 0.891 fm. Thus taking 
the midpoint r p =0.883 fm, we obtain 
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(w), =0-841(6) fm, (65) 

which exactly coincides with the new proton size 1 (i.e. r p =0.S4 184(67) fm). 

4.4. 2Si/2-2P 1 / 2 Lamb shift in He + . 

Modern computed value of this shift is equal to [3] 

Lf s _ 2P =\A 041.46(3) MHz, (66) 

which after the correction (50) becomes 

(fLzP = {r"ifL H 2 l 2P =14 042.21(3) MHz (67) 

(where Y„l=\ .0000266 for He + ). The result of measurement of the Lamb shift by an anisotropy 
quenching method reported in [25], is 

teL 2 , =i4042 - 52 ( i6 ) MHz ' ( 6s ) 

which disagrees with both estimations (66) and (67). 

The obvious discrepancy of the experimental value (68) and QED prediction (66) stimulated 
further experimental research of the 2S-2P Lamb shift in He + . In course of their work the authors of 
ref. [25] redesigned a photon detector system to eliminate a residual polarization sensitivity of the 
photon detectors, which, in authors' opinion, distorted the result of the previous measurement (68). 
Having implemented this improvement, they reported in [26] a new result 

fe) 2j2? =14041.13(17)MHz, (69) 

which again is in disagreement with the alternative predictions (66) and (67). 

Thus, the performance of new high precision experiments on the subject appears to be highly 
required. 

4.5. 1S Lamb shift in hydrogen. 

Having corrected the classic 2S-2P Lamb shift in PBFT, we are now in the position to complete the 
PBFT corrections to the ground state Lamb shift (54). To the accuracy sufficient for further calcula- 
tions, we adopt that the term SL 2S in eq. (54) is completely determined by the corresponding correc- 
tion to the coefficient A in eq. (61a) for 2S-2P Lamb shift. Thus, we put 8L 2S x SA = {v mn 2 —ijA. 
Hence SL 2S =14.08 kHz. Inserting this value into eq. (52), we obtain 

{L h \ s =L IS + 73.0 kHz. (70) 

Our next goal is to estimate the proton charge radius derived from the ground state Lamb 
shift corrected by eq. (70) via the comparison of calculated and experimental data on 15* Lamb shift 
collected in Table 12.3 of ref. [3]. One should point out that major part of experiments for the mea- 
surement of the 15* Lamb shift in hydrogen has been carried out with standard radiofrequency me- 
thod, whose data are rather widely scattered between the values 8 172 798 kHz and 8 172 874 kHz, 
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with the typical measurement error about 30-50 kHz. In these conditions we select the result of 
the mentioned Table [3] 

L ls =8 172 837(22) kHz. (71) 

obtained with Doppler-free two-photon laser spectroscopy [18], which provides more accurate de- 
termination of both 2S-2P and IS Lamb shifts in comparison with radiofrequency method. 

The same work [3] presents the corresponding theoretical values (L IS ) for two different val- 
ues of proton charge radius: 

Z 1S =8 172 663(6) kHz (r p =0.805(l 1) fm) [23], and (72a) 

L 1S =S 172 81 1(14) kHz (forr p =0.862(12) fm) [24]. (72b) 

According to eq. (70), we have to add 73.0 kHz to the values (72a-b). Hence we obtain 
(Zj 15 =8 172 736(14) kHz (r p =0.805(l 1) fm), and (73a) 

{L b \ s =S 172 884(14) kHz (for r p =0.862(12) fm). (73b) 

Using the parameterization (60) with the corrected data (73a-b), we find the coefficients \A b \ s and 

(B h \ s as follows: (A b ) ls =S 111 727.2 kHz, (B h ) ls =\ 557.58 kHz/fin. This allows us to determine 

the proton charge radius, equating to each other the corrected theoretical value and the experimental 
value (71). This coincidence occurs at 

r p =0.844(22) fm. (74) 

This estimation again perfectly agrees with the value of r p determined through the 2S-2P 
Lamb shift in PBFT (65) and with the proton size determined in ref. [5]. 

5. 2S-2P interval in Li-like uranium 

We emphasize that our approach is well applicable not only to light hydrogenic atoms, but also to 
the entire atomic physics, including the case of heavy multi-charged ions. Thus the corrections of 
PBFT we have introduced are relevant for both light and heavy atoms. The basic reason, distinguish- 
ing the cases of light and heavy atoms from mathematical viewpoint, is that for the latter case the 
parameter Za is not small (for example, for uranium it is about 0.67). At the moment PBFT is for- 
mulated in the form of perturbation theory based on the Dirac and Breit equations. Translation of 
PBFT to the QED mathematical language is also possible, but stays outside the scope of the present 
work. Nonetheless, we can still apply the approach of perturbation theory in the derivation of PBFT 
corrections for heavy atoms too, at least while the theoretical and experimental uncertainties in esti- 
mation of energy intervals in such atoms remain comparably large (from 10" 3 to 10" 4 in relative 
units). At the same time, the introduction of PBFT corrections to the effects of finite size of nucleus 
and its polarization requires a separate research. The same remark is relevant with respect to intere- 
lectronic interaction for Li-like multi-charged ions, which should be processed within PBFT, too. 

Thus a consistent analysis of heavy atoms in the framework of PBFT will be done elsewhere. 
In the present contribution we can consider, at least qualitatively, the 2S-2P interval in heavy Li-like 
uranium measured with the experimental uncertainly about 3-10" 4 (280.59(10) eV [27] and 
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280.52(10) eV [28]), which remains the most precise result in the physics of multi-charged heavy 
ions. The current theoretical estimation of this interval is 280.71(10) eV [29] and exceeds the expe- 
rimental value by about 0.07 %. 

Correcting this theoretical value within PBFT, we remind again that our approach does not 
imply any modifications of QED structure. Rather we introduce the PBFT corrections to the input of 
resulting QED expressions by the replacements m^mb„ , U^y n U, and r = r'/(b n y n ) for the bound 

electron (see ref. [7]). In addition, we take into account that such corrections even for heavy atoms 
cannot be large (maximum few percents in relative units for each energy term), and further deal with 
only the dominant contributions into 2S-2P interval in heavy Li-like U: interelectronic interaction 
with the leading term of one-photon exchange (^=368.83 eV [29]), self-energy and vacuum pola- 
rization (^£=-54.98(7) eV [30]), and finite nuclear size correction (0^^33,35(6) eV [31]). 

Interelectronic interaction is described by the Hamiltonian V ee [32], which is added to the 
Coulomb interaction between the nucleus and electron, so that we obtain 

H = ~— + V ee = - — (l - V e/ /Ze 2 ). (75) 
r r 

In PBFT the Hamiltonian of interelectronic interaction has the form of (75) in r' -coordinates. Hence 
one can see that the maximum correction to this interaction induced by PBFT has the value of 

W ee *W ee (l-b n y n )/Z (76) 

due to the scaling transformation (9). For n=2 state of uranium (Z=92), the parameter 
(l - b mn b Mn y mn y Mr )l Z =6.3- 10~ 4 , and the correction to the interelectronic interaction becomes 

SW ee =368.83-6.3- 10" 4 =0.23 eV. (76a) 

The sign of 5W ee is positive, because the transformation (9) signifies that at the semi-classical level 
the increase of the radius of the orbit of IS electrons in Li-like uranium, which, in turn, does in- 
crease the screening Coulomb effect on 25* and 2P states of this atom. 

The correction SWse to self-energy and vacuum polarization within PBFT is determined 
straightforwardly via eq. (50a) for the effective nuclear charge Z'= Z - 2 = 90 . Hence 

5W SE = -54.98((l - (Z'af/4} ' 5 - 1)= -3.23 eV. (77) 

The PBFT correction SWfns to finite nuclear size emerges due to the fact that in the classical 
analogy of PBFT, the electron orbits around the nucleus at larger radius than in the standard theory 
and thus, it is less sensitive to variation of nuclear size than in the standard approach. In order to 
estimate this effect numerically, we use the approximate expression for nuclear size effect presented 
in [31], which in the general form reads: 

A£„ v - t£ nPm = f(n{Z 4mZe 2 R/n 2 f^ , (78) 

where f(n,(Za)) is some function defined semi-numerically, R is the effective nuclear radius, and 
r B = h 2 1 mZe 2 is the Bohr radius. In PBFT, eq. (78) should be modified by the replacement m— >m&„. 

For n=2 state, l/b n =(l-(Za) 2 /4) , and for the numerical value of rhs of eq. (78) -33.35 eV [31], 
we obtain the PBFT correction to the finite nuclear size effect as follows: 
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5W FNS = -33.35 



\\-(ZaflA) 



=2.83 eV. (79) 



Thus the total correction to 2S-2P interval in Li-like U is determined as the sum of eqs. 
(76a), (77) and (79): 



SW^ = SW ee + 5W SE + 5W FNS = -0.17 eV. 

Hence the corrected within PBFT 2S-2P interval in Li-like uranium becomes 

{E h \s-2 P =E 1S - 2P +SW tolal =280.54(10) eV, (80) 

that is in a perfect agreement with both experimental values mentioned above 280.59(10) eV [27] 
and 280.52(10) eV [28]. 

At the same time, at this stage it would be hasty to interpret this result as another achieve- 
ment of PBFT, since the correction to interelectronic interaction has been evaluated only qualitative- 
ly. Rather this result demonstrates a heuristic potential of PBFT in its application to heavy atoms. 
Thus, it seems very interesting to extend consistently the pure bound field theory to the case of 
heavy atoms, especially with the appreciable success of this theory demonstrated above for the light 
hydrogenlike atoms, including positronium. 

6. Conclusion 



In this paper we verify the Pure Bound Field Theory (PBFT) based on re -postulated Dirac and Breit 
equations at the scale of hyperfine contributions to the atomic energy levels. On this way we consis- 
tently considered the fine structure corrections, as well as corrections to hyperfine spin-spin interac- 
tion and corrections to the Lamb shift. 

We have demonstrated that the fine structure corrections have the order of magnitude 

mc 2 {Zaf m/M and scales as n' 5 and n' 6 . Hence they are practically significant only for 15* state of 
hydrogenlike atoms, in particular, in the re-estimation of ground state Lamb shift in the hydrogen 
atom. For IS-2S interval in positronium, there appears an additional component of correction due to 
the appropriate modification of annihilation term in PBFT, and both corrections completely elimi- 
nate the available discrepancy between theoretical value and experimental data for \S-2S interval. 

The corrections to be brought by PBFT to the spin-spin interval occur negligible for the hy- 
drogenlike atoms with m«M (for hydrogen and heavier atoms, eq. (44), it is due to the negligible 

value of the correcting factor {Zaf ^wM W s _ s in comparison with calculation uncertainly; for 

[M + m) 

muonium, eq. (45), it is due to cancellation by the same PBFT correction in the magnetic moment to 
mass ratio value derived from Zeeman effect). At the same time, such PBFT correction acquires 
significant value for the spin-spin interval of 15* state of positronium (eq. (46)), where the correction 
of PBFT removes the discrepancy between theoretical and experimental results. 

The PBFT correction to the Lamb shift emerges due to the replacements (5a-c) in corres- 
ponding QED equations. Such a correction is directly applicable to the classic 2S-2P Lamb shift (eq. 
(50)), where for the hydrogen atom we observe an exact coincidence of corrected by us theoretical 
value of proton charge radius r p =0. 841(6) fm (eq. (65)) with the latest result of measurement via the 
2S-2P Lamb shift in muonic hydrogen r p =0.84 1 84(67) fm [5]. The obtained correction to 2S-2P 
Lamb shift contributes to the corresponding correction to 15* Lamb shift via eq. (54), along with the 
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correction of PBFT to 1S-2S interval (34b). Introducing the PBFT corrections to the 15* Lamb 
shift in hydrogen, we obtained the proton charge radius r p =0. 844(22) fm (eq. (74)), which practical- 
ly coincides with the result yield by the classic 2S-2P Lamb shift. 

We have considered 2S-2P interval in Li-like uranium and shown that PBFT also gives a bet- 
ter coincidence of theoretical and experimental values within the range of their uncertainty. 

In Table 1 we summarize the results of QED without and with the corrections we introduced, 
in comparison with corresponding experimental data. These data completely support our principal 
idea to modify the Dirac equation for non-radiative EM field of bound electron, which, in its turn, 
lead to further modifications of equations of the atomic physics. These modifications induce correc- 
tions into the effects to be not directly related to each other, but characterized by the same final re- 
sult: practical elimination of deviations between theory and experiment. 

We emphasize that the fine structure correction of PBFT in positronium and corrections to 
spin-spin interaction result from the appropriate modification of Breit equation in PBFT. In the 
analysis of radiative corrections to the atomic energy levels, we consider PBFT as a complementary 
to QED, and no modifications of QED core structure are implied. Nonetheless, the corrections of 
PBFT to QED results do emerge due to the replacement of electric potential Uby y mn Um the result- 
ing QED expressions. 

Further, it would be fair to bring up that this work is initiated based on an idea of the third 
author that the rest mass of any object bound to a given field should be decreased as much as the 
mass equivalent of the "static binding energy" coming into play (and this, for classical particles, 
already at rest) [33, 34]. It is worth to note that, via such an approach, he was able to predict the de- 
cay rate retardation effect of bound muons [35]. However, a detailed discussion of this idea falls 
outside the scope of the present paper. 
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Table 1 

Correction of QED results within PBFT in comparison with corresponding experimental (recommended) values for those parameters, where a high measuring precision has 
been achieved 



Parameter 


QED result 


Result corrected by us 


Experimental value(s) 


Ratio 

"deviation/uncertainty" 
before correction 


Ratio 

"deviation/uncertainty" 
after correction 


15-25 interval in 
positronium 


1 233 607 222.2(6) MHz 


1 233 607 215.8(6) MHz 
(sub-section 4.1, eq. (57)) 


1 233 607 216(1) MHz [2] 


3.0- 


<1 


Spin-spin splitting 
in positronium 


203 391.7(8) MHz 


203 386(1) 
(sub-section 4.2, eq. (59)) 


203 389(2) [18] 
203 387(2) [19] 


1.5 

2.5 


-1.5 

<1 


Spin-spin splitting 
in muonium 


4 463 302.88(55) 


remains non-corrected 
(sub-section 3.2, eq. (45)) 


4 463 302.78(5) [20] 


<1 


<1 


Proton charge radius (25- 
2P Lamb shift in hydrogen) 


0.876(6) fin 


0.841(6) fin 
(sub-section 4.3, eq. (65)) 


0.84184(67) fin [5] 


5.7 


<1 


Proton charge radius (15 
Lamb shift in hydrogen) 


0.876(6) fin 


0.844(22) fin 
(sub-section 4.5, eq. (74)) 


0.84184(67) fin [5] 


5.7 


<1 


2S-2P interval in Li-like 
uranium 


280.71(10) eV 


280.54(10) eV 
(section 5, eq. (80)) 


280.59(10) eV [26] 
280.52(10) eV [27] 


1.0 
1.9 


<1 
<1 
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